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ABSTRACT 

We present spectroscopic and imaging observations of the active T dwarf 2MASS 1237+6526, 
intended to investigate the emission mechanism of this cool brown dwarf. The Ha emission hne 
first detected in 1999 July appears to be persistent over 1.6 years, with no significant variation 
from logiQiLfja /^boi) — —4.3, ruling out flaring as a possible source. The relatively high level of 
emission in this object appears to be unique amongst observed late-L and T dwarfs. One of our 
spectra shows an apparent velocity shift in the Ha line, which could be indicative of an accretion 
hot spot in orbit around the brown dwarf; further confirmation of this shift is required. J-band 
monitoring observations fail to detect any significant variability (e.g., eclipsing events) at the 
±0.025 mag level over periods of up to 2.5 hours, and there appears to be no statistical evidence 
of variability for periods of up to 14 hours. These limits constrain the mass of a hypothetical 
interacting secondary to M2 < 20 Mjup for inclinations i > 60°. While our observations do not 
explicitly rule out the binary hypothesis for this object, it does suggest that other mechanisms, 
such as youthful accretion, may be responsible. 

Subject headings: binaries: general — infrared: stars — stars: activity — stars: individual (2MASSI 
J1237392+652615) — stars: low mass, brown dwarfs — stars: variables: general 



1. Introduction 

A previously undetected population of T dwarfs, 
brown dwarfs exhibiting CH4 absorption bands at 
1.6 and 2.2 /xm (Kirkpatrick et al. 1999; Bur- 
gasser et al. 20G2a), have recently been identified 
by various deep optical and near-infrared surveys 
(Strauss et al. 1999; Burgasser et al. 1999, 2000b,c, 
2002a; Cuby et al. 1999; Tsvetanov et al. 2000; 
Leggett et al. 2000; Geballe et al. 2002). One of 
these objects, 2MASSI J1237392+652615 (Bur- 
gasser et al. 1999, hereafter, 2MASS 1237-1-6526), 
identified in the Two Micron All Sky Survey 
(Skrutskie et al. 1997, hereafter 2MASS), was 
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found to have Ha in emission, with a relative 
luminosity logiQ{LHa/^boi) = —4.3 (Burgasser et 
al. 2000a). While many low-mass stars and young 
brown dwarfs are seen to exhibit this feature, likely 
generated by chromospheric magnetic activity, the 
frequency and relative luminosity of Ha emission 
drops rapidly beyond spectral types M7 V (Kirk- 
patrick et al. 2000; Gizis et al. 2000). Indeed, at 
the time of its detection, the emission of 2MASS 
1237-1-6526 was unique amongst objects later than 
type L5 V. Kirkpatrick et al. (2002) have recently 
reported weak Ha emission in the bright T2 V 
SDSSp J125453.90-012247.4 (Leggett et al. 2000, 
hereafter SDSS 1254-0122). 

In order to explain the activity of 2MASS 
1237-^6526, Burgasser et al. (2000a) have hy- 
pothesized that this object may actually be an 
interacting brown dwarf binary system. In this 
scenario, the binary mass ratio and orbital sep- 
aration are such that the lower-mass secondary 
fills its Roche lobe and steadily loses mass to the 
primary. For a 70 yijup^ primary and mass ratio 

We adopt the definitions Mj^p = 1.90x10^° g = 



q = M2/M1 > 0.07 (i.e., M2 > 5 Mj^p), sustained 
mass loss can occur for physical separations a < 6 
Rjup and orbital periods p < 5.5 hours. The close 
orbit required for this scenario currently rules out 
direct imaging of the binary pair (a < 0'.'0002 as- 
suming a distance of 14 pc; Burgasser et al. 1999); 
however, partial eclipsing could be observed in a 
system with an inclination i > 60°. 

We have obtained followup spectroscopic and 
photometric observations of 2MASS 1237+6526 in 
order to test the binary hypothesis, as well as ex- 
amine the possibility that this object was origi- 
nally observed during a prolonged flare. In §2, we 
present red optical (6300-10100 A) spectral data 
obtained over 1.6 yr using the Keck 10m Low Res- 
olution Imaging Spectrograph (Oke et al. 1995, 
hereafter LRIS), and discuss the behavior of the 
Ha emission line over this period. J-band mon- 
itoring observations obtained over two nights us- 
ing the Palomar 60" Near-Infrared Camera (Mur- 
phy et al. 1995, hereafter IRCam) are presented 
in §3. We discuss how the imaging observations 
constrain the binary scenario in §4. Finally, in §5 
we discuss how these follow-up observations con- 
strain various emission mechanisms for this unique 
object. 

2. Spectroscopic Observations 

2MASS 1237+6526 has been observed on seven 

separate occasions in the red optical regime start- 
ing in 1999 July. Early data, obtained on 1999 
July 16-18 (UT) using LRIS, are described in Bur- 
gasser et al. (2000a). Here, we present our latest 
LRIS observations obtained on 2000 March 5 and 
2001 February 20 (UT). Data acquisition and re- 
duction procedures for the first run are described 
below, and are similar to those used for the 2001 
February data as described in Kirkpatrick et al. 
(2002). 

2.1. Data Acquisition and Reduction 

Conditions during our 2000 March 5 (UT) ob- 
servations of 2MASS 1237+6526 were clear with 
0'.'8-0'.'9 seeing. Because of its optical faintness (R 
> 26; Matthews et al. 1996, Burgasser et al. 1999), 
the target was invisible on the guider camera, so 

9.55x10^"' Mq, Rj„p = 7.15x109 cm = 0.10 R© (Cox 
2000; Tholen, Teiffel, & Cox 2000). 



acquisition into a 1" slit was done by blind off- 
set from a nearby optical source using positions 
determined from 2MASS coordinates. Spectral 
data were obtained using the 400 lines mm^^ grat- 
ing blazed at 8500 A, covering the spectral range 
6300 to 10100 A at 9 A resolution. The OG570 
order blocking filter was used to reduce second- 
order light. Two exposures of 1800s each were 
obtained; however, the data from the first expo- 
sure were fainter than expected, so the target was 
repositioned prior to the second exposure. Obser- 
vations of the DZA4 white dwarf WD 1225-079 
(Green, Schmidt, & Liebert 1986; Kilkenny 1986) 
and the Bl V flux standard Hiltner 600 (Hiltner 
1956; Hamuy et al. 1994) were also obtained for 
calibration using the same instrumental configu- 
ration. 

Science data were reduced using standard 
IRAF^ routines. A series of one-second dark ex- 
posures were median-combined and subtracted 
from the science data to remove detector bias. 
Median-combined quartz-lamp flat-field exposures 
(reflected off of the interior dome) were used to 
normalize the response of the detector. Spec- 
tra were then extracted using the APEXTRACT 
routine. Curvature of the dispersion line was de- 
termined from bright standard star exposures and 
applied to the extraction of the 2MASS 1237+6526 
spectrum. Wavelength calibration was done using 
Ne-Ar arc lamp exposures taken immediately af- 
ter the target observations. We then computed 
a telluric correction from the white dwarf spec- 
trum by interpolating over telluric O2 (6850-6900 
A A-band, 7600-7700 A B-band) and H2O (7150- 
7300, 8100-8350, and 9300-9650 A) bands, and we 
applied this correction to our flux standard and 
2MASS 1237+6526 spectral Finally, flux cal- 
ibration was done by correcting our Hiltner 600 
observations to the spectrophotometric data given 
by Hamuy et al. (1994), being careful to interpo- 
late over telluric features and the Balmer Ha line. 
This correction was then applied to the spectrum 

^Imagc Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, 
which arc operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 

^Despite the presence of Ca II HK absorption in this object 
(Kilkenny 1986; Liebert, Wehrse, & Green 1987), the triplet 
lines at 8498, 8542, and 8662 A were not present and hence 
did not aifect the telluric correction. 
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of 2MASS 1237+6526. 

Reduced spectral data from our 2000 March 
and 2001 February observations are shown in Fig- 
ure 1. Full spectra with flux plotted on a loga- 
rithmic scale are shown in the right panel, while 
a close up of the Ha line on a linear scale is 
shown on the left. A combined spectrum from 
our 1999 July 16 (UT) observations is also plot- 
ted for comparison (Burgasser et al. 2000a). The 
2000 March spectrum is roughly 2.6 times fainter 
than our 1999 July observations, likely the result 
of poor targeting onto the slit due to the blind 
offset and proper motion of 2MASS 1237+6526 
(first observed by 2MASS on 1999 March 13 UT). 
The 2001 February spectrum, on the other hand, 
is brighter by a factor of 1.3; atmospheric trans- 
parency on this night was excellent, and the tar- 
get was visible on the facility guide camera, en- 
abling proper slit placement. Features present 
in the spectra arc noted, including the strong 
3(1^1,1^3) and 2(//i,z/;5)+2:/2 H2O bands starting at 
9250 and 9540 A (A uman 1967), a possibly weak 
4(z/i,fc'3) CH4 band centered at 8950 A (Dick & 
Fink 1977), and lines of Cs I at 8521 and 8943 A 
(6s^5i/2^6p^Pi/2_3/2). The dominant feature in 
these spectra is the strongly pressure-broadened K 
I resonance doublet (4s^S'i/2-^4p^Pi/2. 3/2)1 cen- 
tered at 7665 and 7699 A, but extending over 
most of the spectral regime shown. Liebert et al. 
(2000) have shown that the red wing of the K I 
doublet is the primary contributor to the steep 
slope seen from 8000-10000 A in the spectrum of 
SDSSp J162414.37+002915.6 (Strauss et al. 1999), 
and this feature is likely responsible for shaping 
the red optical spectrum of 2MASS 1237+6526 as 
well. 

2.2. Hq Emission in 2MASS 1237+6526 

Ha emission at 6563 A is clearly seen in all of 
the spectra of Figure 1, and its presence in our 
recent data suggests that emission is persistent in 
2MASS 1237+6526 over a period of at least 1.6 
years. We have investigated the temporal behav- 
ior of Ha emission by measuring the line flux in all 
of our spectra to date. Fltix data were integrated 
within 20 A of the emission peak after subtracting 
off the neighboring pseudo-continuum, the scatter 
of which was used to estimate uncertainties. Cor- 
rections have also been computed to compensate 
for slit losses, by scaling the individual spectra to 



match data from 1999 July 16 (UT) in the 8000- 
9500 A region. The initial and corrected line fluxes 
are listed in Table 1. 

Figure 2 plots these line fluxes versus UT time. 
All of the corrected values are in excellent agree- 
ment with the mean line flux from our 1999 July 
observations, (7.0+0.5) x 10""'^'^ erg cm~^ s""'^, with 
the exception of data taken on 1999 July 17 (UT). 
As reported by Burgasser et al. (2000a), these data 
were of extrenicily low signal-to-noise (again, likely 
due to poor slit placement), and it is possible that 
the formal errors are too low. Excepting this data 
point, we find that the Ha emission in 2MASS 
1237+6526 is very stable, firmly ruling out any 
kind of flaring event. The weighted mean value 
over all observations is fna = (7.44+0.08) x 10"-^'^ 
erg cm~^ s~^, consistent with previous results. 

2.3. Comparison to Other Cool Stars and 
Brown Dwarfs 

To gauge how the strength of Ha emission 
in 2MASS 1237+6526 compares to other active, 
cool dwarfs, we computed the relative Ha lu- 
minosity for each spectrum, log2o(L//a/L6o() = 
logio(///a//boi), where 

logio/bo; = -0.4rn.6o; - 4.61 (1) 

in units of erg cm^^ s^^ (Drilling & Landolt 2000), 
and niboi = rrij + BCj is the apparent bolo- 
metric magnitude. Using a revised J-band mag- 
nitude for 2MASS 1237+6526 from the 2MASS 
Second Incremental Release Catalog and BCj = 
2.09+0.10 (Leggett et al. 1999, determined from 
measurements of Gliese 229B^), we derive niboi = 
18.12+0.14 and a mean \ogiQ{LHa/^boi) = —4.3, 
consistent with the results of Burgasser et al. 
(2000a). 

We compare this value to the emission ratios 
of late-M, L, and T dwarfs in Figure 3. Using 
spectral data from Kirkpatrick et al. (1999, 2000, 
2002) and Burgasser et al. (2000a, 2002b), we com- 
puted Ha line fluxes directly from the spectra as 
described above, with bolometric fluxes for late-M 
and L dwarfs computed from 2MASS J-band pho- 
tometry and BCj values from Reid et al. (2001). 

'Gliese 229B and 2MASS 1237+6526 have the same spectral 
type, T6.5 V (Burgasser et al. 2002a); hence, differences in 
the bolometric correction are likely to be negligible. 
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For T dwarfs, wc used 2MASS J-band photome- 
try and BCj = 2.09, yielding an upper limit to 
the true bolometric correction for the early- and 
mid-T dwarfs (Burgasscr et al. 2002a), and hence 
a lower limit to the bolometric flux; i.e., overesti- 
mating the relative Ha luminosity. Detections are 
plotted as open symbols, while upper limits (de- 
termined from pseudocontinuum noise levels) are 
indicated by downward arrows. For comparison, 
wc have included M and L dwarf samples from 
Hawley, Gizis, & Reid (1996, crosses) and Gizis 
et al. (2000, diamonds). The trend of decreas- 
ing activity beyond M5- M7 V reported by Gizis 
et al. (2000) is readily apparent in this plot, and 
appears to continue well into the T dwarf regime. 
Indeed, only two objects have reliable detections 
beyond spectral type L5 V, the T2 V SDSS 1254- 
0122, with logio(L/fa/Lbo;) = -5.7, and 2MASS 
1237-1-6526. The former object is over 25 times 
less luminous in Ha then 2MASS 1237-^6526, and 
upper limits for the other T dwarfs do not exceed 
logio(LHa/L&oi) = -5.2. 

Burgasser et al. (2000a) compares 2MASS 
1237-1-6526 to another unusually active cool dwarf, 
the M9.5 Ve PC 0025-^0047 (Schneider et al. 
1991), which is indicated in Figure 3 by the solid 
square. This object, with logjolLffa/Lhoi) = —3.4, 
is also more than ten times brighter in Ha than 
other objects in its class, although it is not ex- 
tremely active in comparison to earlier emission 
M dwarfs (Hawley, Gizis, & Reid 1996). PC 
0025-1-0047 also shows long-term stability in its 
emission (Schneider et al. 1991; Mould et al. 1994; 
Martin, Basri, & Zapatero Osorio 1999). These 
common characteristics provide further evidence 
that emission processes in these two objects may 
be related. 

2.4. A Possible Velocity Shift in the 2000 
March Spectrum? 

Examination of the spectra of Figure 1 also 
shows an apparent blueshift in the Ha line in our 
2000 March data as compared to the other epochs. 
We measured a line center deviation of 3. 2 ±0.5 A 
in this spectrum with respect to the 2001 February 
data, where the uncertainty is derived from the 
scatter of telluric OH emission line wavelengths 
between these two data sets. This deviation cor- 
responds to a velocity of 150±25 km s~^, which 
cannot be due to differences in the Earth's rota- 



tional and orbital motion between the two epochs, 
amounting to only 2.4 km s~^. 

As with our initial detection of Ha, we carefully 
scrutinized the reality of this subtle (1-2 pixel) line 
shift. Figure 4 shows the Ha region in the raw 
image data from 2000 March and 2001 February. 
Telluric OH emission lines in the 6563 A region ap- 
pear as vertical bands, while the bright spots are 
Ha emission in 2MASS 1237-^6526; note the lack 
of continuum in this region. There is a slight off- 
set in the emission between these two epochs with 
respect to the telluric lines, confirming the spec- 
tral results in Figure 1. However, this shift could 
be due to poor centering of the object on the slit. 
Both individual exposures taken on 2000 March 5 
(UT) show a similar offset, despite repositioning 
of the object. The second exposure was brighter 
and probably better positioned, but was still 2.6 
times fainter than our 1999 July results; hence, 
it too may have been misaligned. On the other 
hand, seeing on the nights of 2000 March 5 and 
2001 February 20 (UT) was roughly equivalent to 
the width of the slit employed, so significant shift- 
ing of the spectrum due to poor slit placement is 
less likely. We were unable to perform an inde- 
pendent check using the weak 8521 and 8943 A 
Cs I lines due to the low signal-to-noise of both 
spectra. We therefore conclude that there; is some; 
evidence for a velocity shift in the Ha line which 
must be confirmed with higher resolution data. 

Assuming its reality, the Ha line shift (which 
appears to be asymmetric toward the blue) could 
arise from a single hot spot on the surface of a 
rapidly rotating brown dwarf. The velocity of the 
line is comparable to but less than the break-up 
speed for a substellar object, roughly 330 (520) km 
for a 30 (75) Mj„p brown dwarf. In this case, 
however, the integration time of each individual 
observation made (typically 20 30 minutes) would 
be a significant fraction of the rotation period (no 
less than about 50 minutes) , and one would expect 
the line to be broadened rather than shifted. A 
more reasonable source would be a hot spot in or- 
bit around the brown dwarf, possibly arising from 
an optically thin accretion disk. Discrete regions 
in an edge-on Keplerian disk would have line-of- 
sight radial velocities as high as 150 km s~^ at 
separations of roughly 5 (12) Rj„p for a 30 (75) 
Mj„p central source. The hot spot would have to 
be small, due to the lack of blue continuum seen in 
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at shorter wavelengths (Burgasser ct al. 2000a). It 
is worth noting that the magnitude and direction 
of the velocity shift is similar to that seen in some 
solar flaring events during an impulsive evapora- 
tion phase (Antonucci et al. 1985); however, we do 
not see the typically brighter redshifted condensa- 
tion component (Fisher, Canfield, & McClymont 
1985). 

3. Imaging Observations 

In addition to our spectroscopic data, we 
searched for photometric variability in 2MASS 
1237+6526, for the purpose of detecting an eclipse 
by a closely-separated companion. 

3.1. Data Acquisition and Reduction 

J-band monitoring observations were obtained 
over two nights, 2000 May 19-20 (UT), using the 
Palomar 60" IRCam. Conditions on both nights 
were generally hazy but clear, with average seeing 
of 1'.'2 to 1"8; some scattered cirrus was present 
during the observations of May 20. Data were ac- 
quired in sets of five 45s exposures dithered by 10" 
on May 19 and 15" on May 20 in a square pattern. 
A total of 155 and 165 images were obtained, span- 
ning periods of 2.50 and 2.57 hours on May 19 and 
20, respectively (Table 2). 

Science data were initially pairwise difFca'enced 
to eliminate bias and sky background, and then di- 
vided by a normalized flat-field image, constructed 
by differencing median- combined sets of lamp-on 
and lamp-off dome observations made on May 19. 
Bad pixels were identified in the flat-field image 
and corrected in the science images by interpo- 
lation. The four separate readouts of the IRCam 
NICMOS-3 HgCdTe array occasionally introduced 
additive quadrant biases that we removed by sub- 
tracting off the median from each quadrant. Fi- 
nally, individual images in each set were shifted 
and coadded to produce a final reduced image 
(Figure 5). 

3.2. Differential Photometry and Analysis 
of the Time Series 

For each image, we computed aperture photom- 
etry (/i) for 2MASS 1237-^6526 and other sources 
in the field. Flux within a 2'.'5 (4.04 pixel) ra- 
dius aperture was integrated after subtracting the 
mean background measured in a 6 to 12" annulus 



around each object. Uncertainty estimates of the 
photometry, Sfi, included contributions by shot 
noise in the source and variations in the back- 
groimd (essentially flat field uncertainties). In- 
strumental magnitudes, rrii = — 2.51ogiQ/j, were 
then computed with corresponding uncertainties, 
CTi = 1.086(-^). The scatter in instrumental mag- 
nitudes for 2MASS 1237-^6526 was roughly 0.03 
mag on May 19 and 0.01 mag on May 20, with 
formal photometric uncertainties of ±0.025 mag. 

To compute a relative magnitude light curve, 
we first examined sources in the field for variabil- 
ity, using the procedure outlined in Bailer- Jones & 
Mundt (2001). We first determined relative mag- 
nitudes and uncertainties for each source: 



- U) (k) ref 

m) = TO ■ — TO, 



(2) 
(3) 



where the reference magnitudes and uncertainties 
were computed as 



K 



TO, 



r^ = -2.51ogio|^/i 



(k) 



ref^ 



K 



(4) 



(5) 



summing over the K other sources in the field. 
After subtracting off the mean relative magnitude, 

we computed the reduced for each source to de- 
termine if photometric variations were significant; 



1 " 



(6) 



here, n is the total number of images. We con- 
strained calibrator sources to have < 0.5, 
corresponding to a probability of variability P 
< 1% for n = 30. Two calibrators present in 
images from both nights passed this constraint, 
2MASSI J1237367+652600 (Calibrator 1) and 
2MASSI J1237372+652536 (Calibrator 2). These 
sources have 2MASS J magnitudes of 15.07±0.05 
and 15.34±0.06, respectively, and are indicated on 
Figure 5. 

Relative magnitudes were then computed for 
2MASS 1237+6526 using these two cahbrator 
stars. The resulting time series is shown for both 
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nights in Figure 6. The mean relative magnitude 
has been subtracted off in both datasets. Time 
series for the differenced caUbrator magnitudes 
{mf^ — mp'') arc also shown (offset) for compari- 
son. Scatter in the relative magnitudes of 2MASS 
1237+6526 was 0.020 and 0.014 mag on May 19 
and 20, respectively, less than the formal photo- 
metric uncertainties. Corresponding values are 
0.505 and 0.314, which are greater than those of 
the calibrator stars (0.271 and 0.237, respectively), 
but nevertheless reject variability (at levels above 
our detection threshold) at the > 99% confidence 
level. 

To examine to possibility of a modulated signal 
(e.g., multiple eclipsing events or geometric mod- 
ulation), we constructed periodograms for each 
night's data using the formalism of Lomb (1976) 
and Scargle (1982) as summarized in Press et al. 
(1992). Given the duration of the observations, 
we arc sensitive to periods ranging from roughly 
0.1 1.3 hr on each night. The computed pe- 
riodograms (oversamplcd by a factor of 10) for 
both 2MASS 12374-6526 (black hues) and cali- 
brator stars (grey lines) are shown in Figure 7. 
May 19 and 20 data are identified by solid and 
dashed lines, respectively. None of the peaks in 
the 2MASS 1237+6526 periodograms are signifi- 
cant, as they do not exceed the 50% false alarm 
probability level and are roughly equal in magni- 
tude to insignificant peaks seen in the calibrator 
periodograms. 

Variability over two nights also appears to be 
absent. There is a 0.024 mag offset in the mean 
relative magnitudes between the two nights; this 
is greater than the difference between the mean 
calibrator values (0.003 mag), but comparable to 
the photometric uncertainties. We compute Xr = 
0.593 over both nights, again rejecting variability 
at the > 99% confidence level. Because of the 
large temporal gap in the time series data, pe- 
riodicity analysis requires more sophisticated al- 
gorithms (e.g., the CLEAN algorithm; Roberts, 
Lehar, & Dreher 1987); however, given the lack of 
statistical evidence for variability, further analysis 
appears to be unwarranted. 

Given the null results from these statistical 
tests, we conclude that there is no short-term vari- 
ability in 2MASS 1237+6526 at the ±0.025 mag 
level over periods of 0.1-2.5 hours and no evidence 
of periodicity for periods of up to 1.3 hours. There 



also appears to be no compelling indication of vari- 
ability at this level for periods of up to ~ 14 hours. 

4. The Binary Hypothesis Revisited 

In order to determine what constraints our null 
variability results impose on the interacting bi- 
nary hypothesis, we have updated the basic the- 
ory described in Burgasser et al. (2000a) by com- 
parison with theoretical models from Burrows et 
al. (1997). Our observables are (1) a period 
limit from the duration of monitoring observa- 
tions, (2) a magnitude limit at which we are sen- 
sitive to variability, and (3) the properties of the 
primary. We therefore seek to derive a relation 
= where p and i are the period 

and inclination of the system, and Mi the mass 
of the primary. For this discussion, we normalize 
masses to Mj„p and distances to Rj„p, and assume 
a circular orbit and negligible angular momentum 
loss (e.g., to gravitational radiation or disk inter- 
action). The period (in hours) is determined from 
the standard Keplerian relation: 



= 47r' 



G{Mi + M2) 



= 8.79 



(7) 



in Jupiter-based units. In the case of an interact- 
ing binary, the secondary fills its Roche lobe, and 
hence its radius R2 is determined by 



R2 = Rl = d- 



0.49g3 



0.6g3 +ln(l + g3) 



(8) 



(Eggleton 1983). Brown dwarf radii, however, 
are fairly constant from 5-70 Mj„p (Burrows & 
Liebert 1993), and are well-described by the rela- 
tion (Zapolsky & Salpeter 1969; Stevenson 1991): 



R = ai 



ttM- 



(I + I.8M-3) 



4 ) 
3 



(9) 



where ai is an age-dependent scale factor (essen- 
tially a correction to partial degeneracy) chosen to 
match the interior models of Burrows et al. (1997). 
Because R2 depends only on the primary mass and 
mass ratio q, both a and p are functions of Mi and 
q exclusively. Therefore, these two parameters de- 
fine the dynamics of the system. 

To determine the observable change in J-band 
magnitude caused by a nearly edge-on eclipse, we 
assumed a simple geometry of overlapping spheres 
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along the lino of sight, as diagrammed in Figure 
8. In general, because dlnR/dlnM « — ^ < 
for brown dwarfs, the lower-mass secondary has 
a larger radius (i?2 > Ri) for M2 > 5 Mj„p. The 
relative area of the primary. A, obscured by the 
secondary during maximal eclipse is: 

7rA = 6li-0.5sin26li + (^)2(6'2-0.5sin6'2) (10) 
where 

iiisin6li = ii2sin(92 (11) 

Ri cos 61 + R2COS62 = a cos i, (12) 

with 9 1 and 02 defined in Figure 8. These equa- 
tions allow us to uniquely determine A in terms 
of Ml, q, and i. Note that no eclipse occurs when 
cosz > (i?i + i?2)/a, while a complete eclipse oc- 
curs for cosi < (-R2 — Ri)/a. 

The observed magnitude difference at maximal 
eclipse is then 

AJ = -2Mog,o{^^^y^^) (13) 

where /3 = F2/F1 is the flux ratio of the sec- 
ondary to primary at J-band. Wc can determine 
the change in J-band magnitude as a function of 
Te// empirically from the two known companion 
T dwarfs, Ghcsc 229B (Nakajima et al. 1995) and 
Gliese 570D (Burgasser et al. 2000b). The former 
object is a T6.5 V with Mj = 15.52±0.06 (Leggett 
et al. 1999) and Te// - 950 K (Allard et al. 1996; 
Marley et al. 1996; Tsuji et al. 1996); the latter ob- 
ject is a T8 V with Mj = 16.47±0.07 (Burgasser et 
al. 2000b) and Te// - 800 K (Geballe et al. 2001). 
These values yield the very approximate relation 
dJ/dT = —0.0063 mag/K, which may only be ac- 
curate over a limited Te// rangc"'^". Interior mod- 
els from Burrows et al. (1997) show a power law 
relation between Te// and mass which is roughly 
dependent on age (r). A fit to the models yields 
the relations: 




88.7M-0-43 r = 0.5 Gyr 
73.4M-0-44 ^ ^ ;^_Q K/Mj, 
45.4M-0-44 T = 5.0 Gyr 



(14) 



This rough relation also ignores the effects of gravity and 
metallicity, which could be substantial. 



Combining Equation 14 with the expression for 
dJ/dT yields the secondary flux ratio: 

-2.5/offio/3 = J^I^AAf (15) 
= a2x(l-g) 

where a2 is a constant dependent on age and 
primary mass. We have assumed that 2MASS 
1237-F6526, a T6.5 V similar to Gliese 229B, has a 
Te// ~ 950 K, and masses 23, 30, and 55 Mj^p for 
T = 0.5, 1.0, and 5.0 Gyr, respectively (Burrows et 
al. 1997). The various parameters are summarized 
in Table 3. 

Using the above equations, wc have c;ak;iilated 
the expected magnitude drop (AJec;) at minimum 
light for an eclipsing, interacting brown dwarf bi- 
nary with Tprimary = 950 K, 0.5 < T < 5.0 Gyr, 
<q<l, and 40° < i < 90°. Results are shown in 
Figure 9 for r = 0.5 (solid line), 1.0 (dot-dashed 
line), and 5.0 (dashed line) Gyr, and (from top to 
bottom) i = 90°, 80°, 70°, and 60°. Only values 
of AJeci for which q < 0.63 and M2 > 5 Mjup 
are displayed, corresponding to the limits for sus- 
tained mass loss (Burgasser et al. 2000a). Signif- 
icant drops in the observed flux are predicted by 
this model, particularly for i > 80°, where fluxes 
at maximum eclipse over 0.7 mag fainter than the 
nominal brightness are possible over the full range 
of periods shown. These flux deviations are due 
primarily to the substantial brightness ratio be- 
tween the secondary and primary, which is partic- 
ularly important at high inclination angles where 
obscuration is most complete. At lower inclination 
angles, the limiting factor is the orbital separation, 
which is larger for lower mass secondaries. No 
eclipse variations are predicted to occur for i < 50° 
over the range of ages examined. Note that sig- 
nificant systematic errors may be present for the 
lowest mass ratios (i.e., low secondary Te//) due 
to the rough estimate of (3 in Equation 15; how- 
ever, our observations generally sample only the 
highest mass ratios due to the period cutoff. 

The predicted AJec; from these models are gen- 
erally much larger than our observational limits 
(indicated by the hatched region in Figure 8). The 
constraints imposed by the non-detection of vari- 
ability are summarized in Table 3 for i > 60°; in 
this case, q > 0.70, 0.53, and 0.37, corresponding 
to secondary masses M2 > 16, 16, and 20 Mj„p for 
r = 0.5, 1.0, and 5.0 Gyr, respectively. Note that 
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for the 0.5 Gyr case, our mass ratio hmit exceeds 
the minimum q required to sustain mass loss. The 
lowest inclination angle eclipses that could have 
been detected by our observations for r = 1.0 and 
5.0 Gyr are imin = 54° and 56°, with correspond- 
ing minimum secondary masses of 20 and 34 M j„p, 
respectively. 

5. Discussion 

Based on the results of the monitoring observa- 
tions, we can confidently rule out an interacting 
binary system for 2MASS 1237-^6526 for i > 60° 
and M2 > 20 Mj„p. We cannot rule out a less 
inclined system, although such a geometry would 
fail to explain the apparent line shift seen in our 
2000 March spectral data. A compromise solu- 
tion would be a moderately inclined system with 
a massive primary and large mass ratio. Conclu- 
sive synthesis of the spectroscopic and photomet- 
ric data, however, requires verification of Ha line 
shift. 

Other emission mechanisms should also be con- 
sidered. 2MASS 1237-1-6526 could be a young, 
and hence very low mass, brown dwarf that is 
still accreting material from a circum(sub)stellar 
disk. If this object were in a weak-lined T-Tauri 
(WLTT) phase, its Ha line could arise from a 
small, optically thin accretion boundary layer, 
consistent with our spectral and photometric ob- 
servations. Assuming a typical WLTT age (30-100 
Myr; Hartmann, Kenyon, & Hartigan 1993) and 
Teff ^ 950 K, the evolutionary models of Bur- 
rows ct al. (1997) predict a mass of 3-12 Mj„p, 
below the Deuterium-burning limit. Note that 
2MASS 1237-F6526 is not associated with any 
known star-forming region, although it could have 
been ejected from its nascent cluster (Reipurth & 
Clarke 2001). The absence of the 1.25 /zm K I 
doublet in this object (I. McLean, priv. comm.) 
marginally supports the possibility of it being a 
low-gravity source, although we cannot rule out a 
purely temperature effect (Burgasser et al. 2002a). 

Our follow-up observations of 2MASS 1237-^6526 
have enabled us to rule out flaring as a viable 
means of emission, and have placed significant 
constraints on the geometry and membership of an 
interacting binary system. Parallax (to determine 
absolute brightness) and space motion measure- 
ments are clearly required to further characterize 



this object. Furthermore, the apparent velocity 
shift described in §2.4, which may prove to be 
a vital clue to the origin of emission in 2MASS 
1237-1-6526, requires confirmation and more de- 
tailed, higher-resolution investigation. The nature 
of this enigmatic brown dwarf may hopefully be 
revealed with further observations. 
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Fig. 1.— Optical spectra of 2MASS 1237+6526 
obtained using LRIS. Observations from 1999 July 
16, 2000 March 5, and 2001 February 20 (UT) 
arc shown, and arc scaled by the correction fac- 
tors described in §2.2. The right panel shows full 
red optical spectra from 6300-10100 A, with flux 
plotted on a logarithmic scale to enhance features; 
spectra are offset by a constant factor. Observed 
bands of H2O and CH4, and lines of Cs I and Ha 
are labelled, along with the pressure-broadened K 
I doublet. The left panel shows a closeup of the 
Ha line on a linear scale, with data offset by an 
additive constant (dotted lines). 

Fig. 2.— Emitted Ha flux from 2MASS 
1237+6526 versus UT time. Circles plot corrected 
line flux values, scaled by matching flux data in 
the 8000-9500 A region to observations made on 

1999 July 16 (UT). The mean emission strength 
fHa = (7.44+0.08) X 10"" erg cm'^ s'^ is indi- 
cated by the dashed line. 

Fig. 3. — Relative Ha luminosity versus spectral 

type for cool dwarfs. The mean relative lumi- 
nosity for 2MASS 1237+6526 is indicated by the 
solid circle, while values calculated from late-M, L, 
and T dwarf spectra are shown as open symbols 
(detected) or arrows (nondetection upper limits). 
Also shown are data for M and L dwarfs from Haw- 
ley, Gizis, & Reid (1996, crosses) and Gizis et al. 
(2000, diamonds), and for PC 0025+0047 (Schnei- 
der et al. 1991, solid square). T spectral types are 
from Burgasser et al. (2002a). 

Fig. 4.— Raw LRIS data of 2MASS 1237+6526 
from 2000 March 5 (top) and 2001 February 20 
(UT, bottom) in the 6563 A region. Vertical 
bands are telluric OH emission lines, while the 
bright spots arise from Ha emission in 2MASS 
1237+6526. 

Fig. 5. — Combination of five, 45 sec J-band im- 
ages of the 2MASS 1237+6526 field obtained on 

2000 May 19 (UT). 2MASS 1237+6526 is indi- 
cated, as are the two comparison stars Cal 1 and 
Gal 2. The field shown is 2'x2', oriented with 
North up and East toward the left. 

Fig. 6.- Relative magnitude time series for 
2MASS 1237+6526 (diamonds) and comparison 
stars (triangles). Data for 2000 May 19 and 20 
(UT) are shown separately. The mean relative 



magnitudes for each night have been subtracted, 
with datasets offset for clarity. 

Fig. 7. — Periodograms of monitoring time series 
data. Black and grey lines trace the periodograms 
of 2MASS 1237+6526 and comparison stars, re- 
spectively, while solid and dashed lines identify 
data from 2000 May 19 and 20 (UT). False alarm 
probability limits of 50% and 10% are indicated. 

Fig. 8. — Geometry of the interacting binary 

model, (a) Orientation along the line-of-sight; a 
is the binary separation, i the inclination. Mi and 
M2 the masses of the primary and secondary, and 
i?i and i?2 their respective radii, (b) Face-on ori- 
entation, showing the relevant angles, 0i and 02, 
used to compute the obscured area of the primary. 

Fig. 9. — Constraints on the binary hypothesis. 
Predicted AJeci during maximal eclipse is plotted 
versus period for q < 0.63; M2 > 5 Mjup] r = 0.5 
(solid line), 1.0 (dot-dashed line), and 5.0 (dashed 
line) Gyr; and (from top to bottom) i = 90°, 80°, 
70°, and 60°. The hatched area covers the pa- 
rameter space probed by the J-band monitoring 
observations. 
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Table 1 
Ha Line Strengths. 







Uncorrected 


Corrected 


Corrected 


UT Date 


tint (s) 


fee" 




log (L/jo/Lfcoi) 


(1) 


(2) 


(3) 


(4) 


(5) 


1999 July 16 


4800 


6.9±0.3 


6.9±0.3 


-4.3 


1999 July 17 


1800 


1.6±0.5 


3.4±1.0 


-4.6 


1999 July IS*" 


1800 


8.4±1.4 


6.8±1.1 


-4.3 


2000 March 5 


1800 


2.9±0.3 


7.4±0.6 


-4.3 


2001 February 20 


2400 


10.6±0.3 


7.9±0.2 


-4.2 



^In units of 10"^''' erg cm"^ 

^Data obtained in the 3800—8600 A regime. 



Table 2 
Log of Imaging Observations. 



UT Date 


UT Time 


Total t (hr) 


No. Images 


Conditions 


(1) 


(2) 


(3) 


(4) 


(5) 


19 May 2000 


05 31 26 - 08 01 55 


2.51 


155 


hazy and clear 


20 May 2000 


06 36 28 - 09 13 25 


2.62 


165 


hazy with some scattered cirrus 



Table 3 

Parameters for the 2MASS 1237+6526 Binary Model. 



T (Gyr) 




Ml (Mj„p) 


a2 (mag)* 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0.5 


1.39 


23 


3.5 


0.70° 


16 


1.0 


1.36 


30 


3.3 


0.53 


16 


5.0 


1.30 


55 


3.1 


0.37 


20 



^Scaling factor determined by comparison to the structure models of Burrows 
et al. (1997). 

^Maximum values permitted by non-detection of variability for -i > 60 ° . 

■^This value is greater than the maximum q — 0.63 required to sustain mass 
loss. 



12 




13 



Ha 



1 n-17 -2 -1- 

1 U erg cm s 



CO 



LA 
cn 



CO 

o 



Q 

Q 



4^ 

cn 
O 
O 
O 

o 



cn 
o 



cn 



CD 
00 

O 




Fig. 2.— 



14 



17 



Relative Magnitude 



o 
b 

o 



o 



o 
o 



o 
b 

Cn 



o 

b 
o 



o 
b 

cn 



O 
O 



O 



o 

M 
O 



O 



O 

O 

o 



cn - 

KJ - 

o r 

hO - 

cn Z 



X 



cn 
O 
cn 



X 



fD 



O 



O 

Q 



b 

cn 
cn 

M 

b 

M 
cn 

b 

cn 



X 



o 



5^ 



— ^ 

— 



Kl 

^ 

— > ^ 



O 
O 

o 



O 



Fig. 6. — 



18 




Fig. 7.— 
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